Not all physical errors can be linear CPTP maps in a correlation space 
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In the framework of quantum computational tensor network, which is a general framework of 
measurement-based quantum computation, the resource many-body state is represented in a tensor- 
network form, and universal quantum computation is performed in a virtual linear space, which is 
called a correlation space, where tensors live. Since any unitary operation, state preparation, and 
the projection measurement in the computational basis can be simulated in a correlation space, it is 
natural to expect that fault-tolerant quantum circuits can also be simulated in a correlation space. 
However, we point out that not all physical errors on physical qudits appear as linear completely- 
positive trace-preserving errors in a correlation space. Since the theories of fault-tolerant quantum 
circuits known so far assume such noises, this means that the simulation of fault-tolerant quantum 
circuits in a correlation space is not so straightforward for general resource states. 



PACS numbers: 03. 67. -a 



Introduction. — Quantum many-body states, which 
have long been central research objects in condensed mat- 
ter physics, statistical physics, and quantum chemistry, 
are now attracting the renewed interest in quantum in- 
formation science as fundamental resources for quantum 
information processing. One of the most celebrated ex- 
amples is one-way quantum computation [3-0]. Once 
the highly-entangled many-body state which is called 
the cluster state is prepared, universal quantum com- 
putation is possible with adaptive local measurements 
on each qubit. Recently, the concept of quantum com- 
putational tensor network (QCTN) which is the 
general framework of measurement-based quantum com- 
putation on quantum many-body states, was proposed. 
This novel framework has enabled us to understand how 
general measurement-based quantum computation is per- 
formed on many other resource states beyond the cluster 
state. The most innovative feature of QCTN is that the 
resource state is represented in a tensor network form [7|- 
9], and universal quantum computation is performed in 
the virtual linear space where tensors live. For example, 
let us consider the one-dimensional open-boundary chain 
of N qudits in the matrix product form 



becomes 



d-l 



d-l 



|*) = ••• E (L\A[k N ]...A[k 1 ]\R)\k If ,...,k 1 ), (1) 
fc 1= o k N =a 

where {|0),...,|rf — 1)} is a certain basis in the d- 
dimensional Hilbert space (2 < d < oo), \L) and \R) are 
D-dimensional complex vectors, and {^[0], ...,A[d — 1]} 
are D x D complex matrices. Let us also define the pro- 
jection measurement M.e,4> ° n a single physical qudit by 
M s ^ = {\ae^),\0 g ^)A2),...,\d - 1)}, where \a g ^) = 



cos 



10) 



e*> sin f |1) 



|0> - 



oi<t> 



COS || 1) 



< 6 < 7r, and < <j) < 27r. If we do the measure- 
ment Me,4> on the first physical qudit of Eq. ((T) and if 
the first physical qudit is projected onto, for example, 
| ae,ct>} as a result of this measurement, the state Eq. (p} 



£ (L\A[k N ]...A[k 2 ] i^}" 9 ' ],, \R)\k N , ...,k 2 ) ® \a e!<t> ), 

where A[a e ^\ = cos § A[0] +e _i0 sin §A[1]. Then, we say 
"the operation \R) — ¥ \R) is implemented in the 

correlation space". In particular, if A[0], A[l], 9, and 
<f> are appropriately chosen in such a way that Ajag^] 
is proportional to a unitary, we can "simulate" the uni- 



A[a B 



\A[a e ,. 



r\R) of the vector \R) in the virtual 



tary evolution 

linear space where A's, \R), and \L) live. This virtual lin- 
ear space is called the correlation space The core 
of QCTN is this "virtual quantum computation" in the 
correlation space. If the correlation space has a sufficient 
structure and if A's, \L), and \R) are appropriately cho- 
sen, we can "simulate" universal quantum circuit in the 
correlation space 0-H, [Io|, [H[ ■ 

For the realization of a scalable quantum computer, a 
theory of fault-tolerant (FT) quantum computation 12- 
[l6j is necessary. In fact, several researches have been 
performed on FT quantum computation in the one-way 
model 0, [l7l - [2H . However, there has been no result 
about a theory of FT quantum computation on general 
QCTN [22|. In particular, there is severe lack of knowl- 
edge about FT quantum computation on resource states 
with d > 3. It is necessary to consider resource states 
with d > 3 if we want to enjoy the cooling prepara- 
tion of a resource state and the energy-gap protection of 
measurement-based quantum computation with a physi- 
cally natural Hamiltonian, since no genuinely entangled 
qubit state can be the unique ground state of a two-body 
frustration- free Hamiltonian [26j j . 

One straightforward way of implementing FT 
quantum computation on QCTN is to encode phys- 
ical qudits with a quantum error correcting code: 
Y^t'lo ■■■Y^t~io( L \ A l k N}---^k 1 }\R}\k N , ...^h), where 



2 



|fcj) (i = 1, ...,N) is the encoded version of |fc,) (such as 
1 6) = |000) and |1) = Jill), etc.) In fact, this strategy 



Refs. 



21[ for the one-way model 



was taken in 
(d = 2), and it was shown there that a FT construction 
of the encoded cluster state is possible. For d > 3, 
however, such a strategy is difficult, since theories of 
quantum error correcting codes and FT preparations 
of the encoded resource state are less developed for 
d > 3. Furthermore, if we encode physical qudits with a 
quantum error correcting code, the parent Hamiltonian 
should no longer be two-body interacting one. 

The other way of implementing FT quantum compu- 
tation on QCTN is to simulate FT quantum circuits in 
the correlation space. Since any unitary operation, state 
preparation, and the projective measurement in the com- 
putational basis can be simulated in a correlation space 
(for a precise discussion about the possibility of the mea- 



surement, see Ref. [111), it is natural to expect that 



FT quantum circuits can also be simulated in a corre- 
lation space. An advantage of this strategy is that the- 
ories of FT quantum circuits for qubit systems are well 
developed _|l 2^ 1 6 1 . In fact, this strategy was taken in 



Refs. @, [HHirfor the one-way model {d — 2). They 
simulated FT quantum circuits on the one-way model. 

In this paper, however, we point out that it is not 
so straightforward to simulate FT quantum circuits in 
a correlation space in general. We first consider the 
simulation of quantum circuits in the correlation space 
of pure matrix product states (post-measurement condi- 
tional states). We show that if d > 3 not all physical 
errors on physical qudits appear as linear completely- 
positive trace-preserving (CPTP) errors in the correla- 
tion space. Since all theories of FT quantum circuits 
known so far assume such noises [l2l - [l6| , this means that 
it is not so straightforward to apply these FT theories 
to quantum circuits simulated in the correlation space of 
pure matrix product states. 

We therefore next consider another way of simulating 
quantum circuits in the correlation space by mixing mea- 
surement results. For the cluster state and the tricluster 
state [22! , such a mixing strategy well works: all CPTP 
errors on a physical qubit (or qudit) can be CPTP errors 
in the correlation space. However, this is not the case for 
other general resource states. As an example, we con- 
sider the one-dimensional AKLT state 2|| 2^|, and see 
that not all physical errors on a physical qutrit can be 
linear CPTP errors in the correlation space even if we 
mix measurement results. This suggests that even if we 
mix measurement results, like the cluster model, the sim- 
ulation of FT quantum circuits in the correlation space 
of a general resource state is not so straightforward. 

Simulation with pure states. — First, let us we con- 
sider the simulation of quantum circuits in the correla- 
tion space of pure matrix product states. We show that 
if d > 3 not all physical errors on physical qudits appear 
as linear CPTP errors in the correlation space. 



Since the MPS, Eq. ([T]), is a resource state for 
measurement-based quantum computation, we can as- 
sume without loss of generality that A[ag^}, A[/3g t< j,], 
A[2], A[3], A[d — 1] are unitary up to constants: 
A[a 9 ^\ = c a U a , A\p e ,$] = cpUp, A[2] = c 2 U 2 , A [3] = 
C3I/3, ... A[d- 1] = c d -iU d -i, where c a , cp, c 2 , ... c d -i 
are real positive numbers, U a , Up, U2, U d ~\ are uni- 



tary operators, and A\flg 



sm 



\A[0] 



-i<t> . 



\A[l] 



This means that any operation implemented in the corre- 
lation space by the measurement Aig,<p on a single phys- 
ical qudit of Eq. (JT]) is unitary. Note that this assump- 
tion is reasonable, since otherwise Eq. (TTJ does not seem 
to be useful as a resource for measurement-based quan- 
tum computation. In fact, all known resource states so 
far [3-i, m m m 31], including the cluster state and 
the AKLT state, satisfy this assumption by appropriately 
rotating each local physical basis. Furthermore, we can 



-, c d ~ 
i,...,fci\ 



such that C = 



,p<z-i 2 _ 



(L\A[k N ]...A[k 1 ]\R)\k N ,...,k 1 ) = 
A ^ \R)\kN, ■■; ki) and we can 



take c a ,cp, c 2 ,. 
1, since J2k 

redefine A[ki\/y/C -t A[ki\. 

Theorem: If d > 3, there exists a single-qudit CPTP 
error £ which has the following property: assume that 
£ is applied on a single physical qudit of Eq. (JI]) . If the 
measurement Aig^ is performed on that affected qudit, 
a non-TP operation is implemented in the correlation 
space. 

Proof: In order to show Theorem, let us assume that 



There is no such £ . 



(2) 



We will see that this assumption leads to the contradic- 
tion that d < 2. 

First, let us consider the state 



(3) 



where U a ++b = | ct> <&j + \b){a\ + I - \a)(a\ - \b)(b\ is the 
unitary error which exchanges \a) and |6), and / is the 
identity operator on a sing le qudit. In Eq. (|3j) , the error 
Ui^2 is applied on the first physical qudit of |\f f ). If wc 
do the measurement Aig^ on the first physical qudit of 
Eq. ([3]), and if the measurement result is |2), Eq. (O 
becomes 



ho , ■ ■ ■ , 



(L\A{k 



N 



A[l] 
\A[l]\\ 



\R)\k N> ...,k 2 )®\2). (4) 



In other words, the operation \R) 



A[l] 



ppjfllfl) is im- 
plemented in the correlation space. By the assumption 
Eq. ([2]), this operation should work as a TP operation in 
the correlation space. Therefore, 



At[l] A[l] 



\\A[l]\\ \\A[1„ 
By taking r\ = ||A[1]|| 2 , we obtain 
At[l]A[l] = r)I 



= I. 



(5) 



(6) 
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Second, let us consider the measurement Me, 4, on 
the first physical qudit of (I® N ^ X ® U ^ 2 V S )\^), 



where s G {0, 1, d - 1}, V = ELo e"^ |p)(p| 



is a unitary phase error, and u> = 2ir/d. If 
the measurement result is \ae,<f>), (e~ 2 "'" cos |j4[2] + 
e -z(0+sw) gm |_,4[i])/^y i s implemented in the cor- 
relation space, where ^ = ||e _2ls " cos §A[2] + 
e -iW+^) s i n |yi[l]||. By the assumption Eq. $2$, this 
should work as a TP operation in the correlation space. 
Therefore, 1 I = cos 2 § A+[2]A[2] + sin 2 § + 
isin6»(e- l W- sw )^t[2]A[l] + e 1 ^ 3 ^ A^[1]A[2]). By the 
assumption that all A's are proportional to unitaries, 
ylt^JAp] = £/, where £ = ||A[2]|| 2 . Furthermore, as 
we have shown, A^[1]A[1] = rjl (Eq. ©). Therefore, 

7 7 = e- 1 ^- 8 ^ A^[2]A[l] + e 1 ^- 8 ^ A^[l]A[2], (7) 

where 7' = ^(7 - £cos 2 § - 77 sin 2 §). 

Finally, let us consider the measurement M.q^ on the 
first physical qudit of (J®^" 1 ® IWEWOI*), where 
i G {0, 1, d — 1}. If the measurement result is \ae,c/,), 
(e~ i4w cosf A[l}+< 



sin I A [2])/ VS is implemented 



in the correlation space, where V~5 = ||e lt "cos|A[l] + 
e -*<*-2ttuj s } n |4[2]||. By the assumption Eq. fltjl, this 
should also work as a TP operation in the correlation 
space. Therefore, 

S'l = e iW+t " ) A t [2]A[l] +e- 4 ^ +tw ^ t [l]A[2], (8) 

where 5' 



sin V 

From Eqs. © 



£sin 2 f -Tycos 2 !). 



and ©, 



3—21(0—80;) 



e/ = [e 

_ e 2i(0+t^) ^ 0. Then, 

2i(4>— suj) _ 2i(4>+tcj)\ 



e 2i ^ +t ^]A^[2]A[l], where e = e -K0-^) 7 

Let us assume that e - 2 *(0- S( ") 
e'l = At[2]A[l], where e' = e/(e- 2i ^- sw -' - e 2 *M>+M). if 
e' = 0, At[2]A[l] = 0, which means A[l] = since A[2] is 
unitary up to a constant. Therefore, e' 7^ 0. In this case, 
A[l] = e"A[2] for certain e" ^ 0, since A[2] is unitary 
up to a constant 
This means 



Hence e" 2 ^-^ - e 2l ^ +t ^ = 0. 



2<f> + {t — s)uj = r s ,t7r, 



(9) 



where r Sj t € {0, 1, 2, 3, ...}. Let us take t = s = 0. Then, 
Eq. (|9]) gives = r ,of (70,0 € {0, 1,2, ...}). Let us take 
s = 1, t = 0. Then, Eq. © gives = § + n, f (ri, & 
{0, 1,2,...}). In order to satisfy these two equations at 
the same time, there must exist ro.o and ri,o such that 
r o.of = 7 + r l,of • If r o,o = n,0> th en = 1/d which 
means d = 00. Therefore r 7^ Then we have 

d = 2/(ro.o — rL,o) < 2, which is the contradiction. ■ 

One might think that if we rewrite the post- 
measurement state Eq. (|4]) as 



(L\A[k N ]...A[k 2 



A[l] 



\\A[l]\R)\ 



\R)\k N ,...,K 2 )®\2) 



and redefine the operation implemented in the correlation 



space as \R) 



the TP-ness is recovered in 



the correlation space. However, in this case, the non- 
lineally appears unless A^[l]/l[l] oc /, and therefore if we 
require the linearity in the correlation space, we obtain 
the same contradiction. 

In short, if d > 3 not all physical errors on physical 
qudits appear as linear CPTP errors in the correlation 
space of pure matrix product states 



Simulation by mixing measurement results. — We have 
seen that if we simulate quantum circuits in the correla- 
tion space of pure matrix product states, not all physical 
errors on a physical qudit can be linear CPTP errors 
in the correlation space. Therefore we must simulate 
quantum circuits in the correlation space with another 
method: we consider the simulation by mixing measure- 
ment results. 

Before studying a concrete example, let us consider the 
effect of a CPTP error on general resource states. Let 
us assume that a CPTP error p — > J2J=i FjpFj , where 

J2J=i FjFj = I> occurs on the first physical qudit of |^). 
If we measure the first physical qudit in a certain basis 
{|rn s )}, we obtain V V ; u ;/•. ,.. /,' ) 2 ® \m s )(m s \, 
where Ej >s = J2k A[k](m s \Fj\k) and 



J2 {L\A[k N ]...A[k r 



W(\i>)) 

(^[k' r }...A^[k' N ]\L)\k r) ...,k N )(k! r ,... ) k' N 



If we trace out \m s ), we obtain ^ s J^j W(Ej, s \R)) 2 . This 
means that the map 



\R){R\^Y, E i>>\ R H R \ E l> 



(10) 



is implemented in the correlation space. Since we can 



show 3J 



. E jt s E h E J,s = '■ the ma P E q- (LTHD is cptp. 

In general, we must do feed-forwarding before tracing 
out |m s ) in order to deterministically implement quan- 
tum gate s in the correlation space. As is shown in 
Refs. |18|, l34( , all CPTP errors on a physical qubit (or 



qudit) can be CPTP errors in the correlation space of 
the cluster state and the tricluster state even if we do 
the feed-forwarding in this mixing strategy. However, it 
is not the case for other general resource states j3f|. As 
an example, let us consider the one-dimensional AKLT 
state, where d = 3, A[0] = -^A, A[l] = iJZ, 

A[2] = ±Z. 

If a CPTP error occurs on the first physical qutrit, and 
if we do the usual measurement-based quantum compu- 
tation on the AKLT state, we obtain [34| 

EEE( 

j p=0 q=0 

^2 W(Q r (si, s r )...Q 2 (si,s 2 )E jiSl \R)) r +l 



(ai,...,a r )eSj 



+h(p, q, r) W(Z r ~ x Ej, 2 \R))r+i ® V(P, q), 



(11) 
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where r is the number of measurements (since the gate 
operation is non-deterministic in AKLT model, we must 
repeat measurements until we near-deterministically im- 
plement the desired gate operation), and Sp q is the set of 
measurement outcomes (si,...,s r ) £ {0, l,2} r \ (2,..., 2) 
such that /(si, s r ) = p and g(si, s r ) = q. Here, 
/(si,...,s r ) = ®i =1 (S ait0 ® 5 Sit i) and g(si, s r ) = 
®Ui(S Si ,i ® S Si}2 ). Also, Q k ( Sl ,...,s k ^,0) = Xe lZe / 2 , 
Q fc (ai,...,fl fc _i,l) - XZe ize ' 2 , Q fc (a ls ...,« fc _ ls 2) = Z 
for si = ... = Sfc_i = 2, and Q k (si, s fc _i, 0) = A, 
Qfc(si,...,s fc -i,l) = AZ, <2fc(si, ...,s fe _i,2) = Z for 
other s^s. We also define h(p,q,r) — Op.o«"g,o if r is 
even, and h(p,q,r) — S Pt o6 qi i if r is odd. Finally, 
77(0, 0), 77(0, 1), 7/(1, 0), and 77(1, 1) are mutually orthogo- 
nal states, which record Pauli byproducts. The first term 
of Eq. (jlip corresponds to the mixture of successful mea- 
surement results (desired rotation is implemented) and 
the second term corresponds to the failed measurement 
results (the desired rotation is not implemented.) 
Equation means that for fixed p and q, the map 

\R)(R\^ 

£ Y,Q(s 1 ,...,s r ,j)\R)(R\QHs 1 ,...,s r ,j) 



+h(p, q, r) Z r - X E j>2 \R) (R^Z^ 1 



(12) 



is implemented in the correlation space, where 

Q(si, ...,S r ,j) = Qr{s\, S r )...Q2(si, S2)Ej tSl . 

For example, let us consider the error with w — 1 and 



F\ — Umo 



10} 



1 



V2 



|Q> |i> 

V2 



HI 



|2)(2| 



where Ums.4, = |ae,^)(0| + IAj^X 1 
if p = 1 and q — 0, we can show that [34( 



r.tp>\k)(k\. Then, 



(si 



-,»r)£ 



, 3 



Q f (si 



,j)Q(si 



*,3) 



if r is even, 
This 



is al + ||1)(1| if r is odd, and f3I + 
where a and /? are certain positive numbers |34 
means that the map Eq. (fT"2|) is not linear CPTP. 

Conclusion. — In this paper, we have studied how phys- 
ical errors on a physical qudit appear in the correlation 
space. We have shown that if d > 3 not all physical er- 
rors can be linear CPTP errors in the correlation space 
of pure matrix product states. We have also shown that 
even if we mix the measurement results, not all physical 
errors are linear CPTP errors in the correlation space of 
general resource states. These results suggest that the 
application of the theories of fault-tolerant quantum cir- 
cuits to the correlation space is not so straightforward. 
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